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A Comments on the proof of Proposition 8.3 



In the proof of Proposition 8.3 the assertion was made that each ^(fQ) was contained in 
some open (relative to A) set Gi such that (83) holds. The case s = d follows from the 
hypothesis that A is contained in a <i-dimensional C 1 manifold since in this case the sets 
Ki may be chosen to be almost clopen relative to A. The case s > d is more delicate and 
requires an argument similar to that used in the proof of Theorem 2.1 given in Section 7. In 
this addendum we provide a sketch of these arguments for the interested reader. 

In this note we let B d (a,r) denote the open ball in R d with center a and radius r, B d (a,r) 
the closure of B d (a,r), and C d (a,r) the boundary of B d (a,r). 



A.l Case 1: s = d 



Let M C R d ' be a C 1 manifold (without boundary), i.e. for each point p G M there is 
some open set (relative to M) V C M containing p, some open ball Bd(q,p), p > 0, and a 
diffeomorphic mapping <p p from Bd(q, p) onto V (i.e. <j) p is a homeomorphism of Bd(q, p) onto 
V and <f)p is continuously differentiable on Bd(q,p)). By composing <p p with an appropriate 
affine mapping on R d we may assume without loss of generality that <f) p (0) = p and that 
<f)' p (0) T 4>' p (0) = I so that </>p(0) is an isometric mapping from onto its range in R d . For 
ijG Bd(0,p), let [x, y] denote the directed line segment from x to y. Then 

My) - M%) = [ U z ) dz = Wiv ~ x ) + [ W z (t)) - W)Mt)- i 1 ) 
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Let e > 0. From the continuity of 4>' p and (JTJ) it follows that there is some < v p < p such that 
4> p is bi-Lipschitz on B d (0, v p ) with constant (1 + e). Since <p p (B d (p, u p )) is a neighborhood 
of p, there is some 5 P > such that B d i(p,5 p ) fl M C (p p (B d (p, u p )). Furthermore, as we 
remark in the proof of Proposition 8.3, for all but a countable collection of < 5 < 5 P , the 
set B d > (p, 5) fl M is T^-almost clopen. 

Now suppose A is a compact subset of M and consider the collection 

C ■= {B d ,(a,5) n A | a e A, 5 < 5 a ,H d (C d >(a, 8) C\ A) = 0}. 

Applying the Besicovitch Covering Theorem we may then choose (ipi,Ki) in the proof of 
Proposition 8.3 where each Ki is of the form (p~ l (B d i(a, 5) fl A) with ipi = (fi a . and ipi(Ki) 
is almost clopen with respect to Ti d \A- We may assume without loss of generality that each 
H d (ipi(Ki)) > for all i. 

As described in the proof of Proposition 8.3 (see the discussion following (82)) we can find 
for 7 > an open and almost clopen set Gi that contains ipi(Ki) and is such that 7i d {Gi \ 
■ipi(Ki)) < 7. Since d{Gi \ ipi(Ki)) C dGi U dipi(Ki) is a set of H^-measure 0, we have 
'HdiGi \ ipi{Ki)) = 7i d {Gi \ ipi(Ki)). By the hypotheses of Proposition 8.3, we may choose 7 
small enough so that 

g a4 {Gi \ UKi)y dls = 9s4(G?y^MT d/s < e/2\ (2) 

By Lemma 3.2 we have 

gUGi) > {9s,mK t )- d ' s +g s4 {G i \i, t {K t ))- d l s ) 

which combined with Q gives (83). 

A. 2 Case 2: s > d 

A closed set A C M d ' is said to be c/-regular if there are positive constants Co, c\ and R such 
that cor d < H d (A fl B d >(a,r)) < c\r d for all a G A and < r < R. The assumption that 
A is a compact subset of B := (Jfc=i 4>k(Gk) where each 0^ is bi-Lipschitz on an open set 
Gk C M. d shows that A C A C B for some rf-regular set A C M d . We may assume without 
loss of generality that g (^(iQ)) > 0, i = 1, 2, Applying the following lemma[T]to each 

4>i(Ki) then shows the existence of some open (relative to A) set Gi D ipi(Ki) such that (83) 
holds. 

Lemma 1. Suppose K C M d ' zs a compact subset of a d-regular set A C M d ' . Assume further 
that s > d and < g sd (K) < g s d (K) < 00. T/ien /or eac/i e > 0, there is some 5 > sitc/j 

for any open (relative to A) set G D K with H d (G \K) < 8. 
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Proof. Note that the topology we use here is the A-relative topology. For convenience, we 
set B(a,r) := S^(a, r) H A for a G A and r > 0. 

Step 1. Let < 5 < (l/2) 4d (we will choose 5 sufficiently small at the end of this proof). 
Suppose G is an open set in A such that K C G and Hd(G \ K) < 5. Since TCd\^ is a 
Radon measure, there is an open set V C A, with G C V such that T-id(V \ G) < 5. Let 
7 := dist(G, A \ V) and note that 7 > (we take 7 = 00 if V = A). 

Step 2. Let {cujv} denote a sequence of optimal s-energy configurations for G and let 
n(u N ) := min Xj€uJN \ {Xi} \x t - Xj\. For p > 0, let 

:= {x; G u;tv I ri(u N ) > pA^ 1/d } 

and let uj]^ := \ oj x n . Then 

£>a0 > ^ (pN- l / d y s = \Zjlj\p- s N s l d (N > 2). (4) 

Let > g s d (K). Then fc > g s ,d{G), and so there is some N\ such that E s (ujn) < kN 1+s ^ d 
for N > Ni and we deduce from (0J) that 

l^i < A:p s (iV > JVO. (5) 



Step 3. Let < f < 1/2 and set 



:= {x, G I 5(x;, z/pA^ 1/d ) n K £ 0} 



and let u;^ := tu^ \ uj 2 n . Let > iVj. be large enough so that pN l l d < min(7, R). For 
N > N 2 and points x i: Xj G uj 1 n with Xj 7^ Xj we have 

fi(x,, pN- 1/d ) C 7 and vpN- 1/d ) n Bfo, z/pA^- 1/d ) = 0. 

Notice that the balls B(xi, upN~^ d )^ x, G Sjy, are pairwise disjoint and contained in V \ K 
and hence we have from the ci-regularity assumption, 

\Z 2 N \c v d p d N- 1 < ^2 7-Ld(B(xi,vpN~ 1/d )) =H d ( (J 5(x i} z/piV- 1 ^) 



< \ K) < 25, 



and so 



Thus 



AT - c (z/p) d ' 
N-luU - _ , . 25 



AM JV l^iVl l^iV 



iV A" — c (z/p) 



> 1 " V " T77T7v7 - (7) 



Step 4. For x, G cj^, there must be some jji £ K such that |xj — y$| < z/pA" Let 
<^jv := {l/i I £ ^at}- For A" > A" 2 , and Xj, Xj G cu^, Xj 7^ Xj, we have 



%i — Vi\ < Z^piV < Z/|Xj — X 



3I> 



and so 

= lUi-Xi + Xi-Xj + Xj-yjl > \xi-Xj\-\yi-Xi\-\yj-xj\ > (l-2u)\xi-Xj\, (8) 



which shows that \uj§\ = \u>^\ 



Step 5. From (jSJ), we have |xj — Xj\ < (1 — 2zv) — yj\ for Xj, Xj G uj 2 n and so, for A" > N 2 , 
we have 



>(i-2-r^f>(i-2^ 



jyl+s/d - V J N l +s /d - V \ N J \u%\ 1+s / d ' 

Taking the limit inferior as A" — > 00 then gives 

g (G) = liminf > (1 - 2z/) s (l - kp s - -™f+'l A g (Jjf) 

where x + = max(x, 0). Choosing p — v — S 1 /^ = gives 

g (G) > (1 - 25 1 /( 4d )) s (l - - 25 1 /7c )+ +sA 2 (K). 



Finally, if e > 0, then we can choose 5 > sufficiently small so that Q holds. □ 
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Abstract 

We investigate the energy of arrangements of ./V points on a rectifiable d-dimensional 
manifold A C M. d that interact through the power law (Riesz) potential V = l/r s , 
where s > and r is Euclidean distance in M. d . With £ S (A, N) denoting the minimal 
energy for such iV-point configurations, we determine the asymptotic behavior (as 
N — > oo) of £ S (A, N) for each fixed s > d. Moreover, if A has positive d-dimensional 
Hausdorff measure, we show that iV-point configurations on A that minimize the 
s-energy are asymptotically uniformly distributed with respect to d-dimensional 
Hausdorff measure on A when s > d. Even for the unit sphere S d C these 
results are new. 

Key words: Minimal discrete Riesz energy, Best-packing, Hausdorff measure, 
Rectifiable manifolds, Uniform distribution of points on a sphere, Power law 
potential 
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1 Introduction 



Determining iV points on the unit sphere S d in M. d+l that are in some sense 
uniformly distributed over its surface is a classical problem that has applica- 
tions to such diverse fields as crystallography, electrostatics, viral morphology, 
molecular modeling, and global positioning. Various criteria (appropriate to 
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the application) for the generation of such points include best-packing, mini- 
mization of energy (e.g., Coulomb potentials), spherical t-designs (cubature), 
maximization of volume of convex polyhedra with N vertices on S d , etc. 

A motivation for the present paper is the analysis of the asymptotic behav- 
ior (as N — > oo) of optimal (and near optimal) iV-point configurations that 
minimize the Riesz s-energy 



over all iV-point subsets {x±, . . . , xn} of S d , where s > is a fixed parameter 
and | • | denotes the Euclidean norm in R d+1 . We remark that as s — > oo, with 
N fixed, the s-energy (1) is increasingly dominated by the term(s) involving 
the smallest of pairwise distances and, in this sense, leads to the best-packing 
problem on S d (cf. [3], [4]). We further note that for s = 1 and d = 2, the 
minimization of (1) is the classic Thomson problem (see e.g. [1], [2], [12], [17]). 

In this paper we investigate the case when s is fixed, s > d, and iV — > oo. 
Significantly our results apply not only to the sphere, but to a class of rec- 
tifiable d- dimensional manifolds embedded in M. d . For such manifolds we de- 
termine, for s > d, the asymptotic behavior of the minimum Riesz s-energy 
as well as the asymptotic distribution of optimal and near optimal iV-point 
configurations. Indeed we shall prove that the latter is given by (i-dimensional 
Hausdorff measure on the manifold and that the minimum iV-point Riesz s- 
energy over the manifold is asymptotically given by C s N 1+s / d when s > d 
and by C^N 2 log iV when s — d. The essential feature of these results (see 
Theorems 2.1, 2.2, and 2.4) is not merely the order of growth of the minimum 
energy as N — > oo, but rather the more delicate verification of the existence 
of the positive constants C s for s > d; a fact which is new even for the case of 
the sphere S d when s > d. Somewhat surprising is the fact that we can give 
an explicit formula for Cd (i.e., for the case s — d) in terms of its Hausdorff 
measure for any compact subset of a d- dimensional C 1 -manifold in W 1 ' (see 
Theorem 2.4 and equation (8)). 

We remark that for < s < d, standard potential theoretic arguments can 
be used for the analysis of the minimum energy points (cf. [9]). However, 
for s > d such methods do not apply. Instead we exploit the scaling and 
translation properties of the energy function together with self-similarity and 
convexity arguments. 

For the remainder of this section we introduce some needed notation and, by 
way of further background, we mention known related results for the sphere 
S d . We devote the next section to the statement of our main results. 

Throughout this paper, ojn = {x±, . . . ,xn} denotes a set of iV (possibly 0) 
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distinct points in W 1 ' . For each real s > the s-energy of lon is given by 

where, as above, | • | denotes Euclidean distance in M. d ' . For A C IR d ' we define 
the iV-point minimal s-energy over A by 

£ S (A,N) := inf £ s (c^). (3) 

By convention, the sum over an empty set of indices is taken to be zero and 
the infimum over an empty set is oo. Hence, S S (A, N) = oo if N is greater than 
the cardinality of A and E s (un) — if N — 0, 1. It is clear that £ S (A, N) = 
£ S (A, TV), where A denotes the closure of A and, furthermore, that £ S (A, N) = 
if A is unbounded. Hence, without loss of generality, we may restrict ourselves 
to the case that A is compact. 

For the unit sphere S d C R d+1 , the asymptotic behavior (as N — > oo) of 
£ s (S d 1 N) is quite different for the three cases (i) < s < d; (ii) s = d; and 
(Hi) s > d. The reason for this is that in case (i), the energy integral 

I.(jm):=[[ —L—dp(x)dii(v) (4) 
JJs d xs d \x — y\ s 

taken over all probability measures \i supported on S d is minimal for nor- 
malized Lebesgue measure Hd(-)\s d /7~Ld(S d ) on S d . However, for s > d, we 
have I s (fJ>) = +oo for all such measures \x. Roughly speaking, as the parame- 
ter s increases, there is a transition from the domination of global effects to 
the domination of more local (near-neighbors) influences, and this transition 
occurs precisely when s = d. 

The following results are known for the above mentioned cases. In case (i), 
classical potential theory yields (cf. [9]): 

Theorem 1.1 If0<s< d, 

£ s (S d ,N) _ ( H d (.)\ s A 

A n 2 - Is { n d (s d ) ) ' (5j 

where I s is defined in (4)- Moreover, any sequence of optimal s-energy config- 
urations (u)^)T c S d i s asymptotically uniformly distributed in the sense that 
for the weak-star topology of measures, 

1 K 7~Ld(-)\s d nr / a \ 



N 



where 5 X denotes the unit point mass at x. 
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For case (it), we have from the results of Kuijlaars and Saff [8] and Gotz and 
Saff [6] the following: 

Theorem 1.2 Let B d := B(0, 1) denote closed unit ball in W 1 . For s = d, 



N 



£ d (s d ,N) _ n d {B d ) _i r(f, 

■>*ooN*\ogN H d (S d ) dv^T(i)' 



and any sequence (ou^) C S" 1 of optimal d-energy configurations satisfies (6). 

(The reader is cautioned that the definition of energy used here differs by a 
factor of 2 from that in [8].) 

Until now, results for s > d have been less complete, describing only the order 
of growth of £ s (S d , N). The following is proved in [8]. 

Theorem 1.3 For s > d, there exist positive constants c\ = ci(s,d), ci = 
c 2 (s, d) such that 

CiN i+s/d < g ^ S d^ N ^ < C2 N 1+a / d , N>2. 



Natural questions that therefore arise for the case s > d are: 

(a) Does the limit 

lim ^ exist? 

(b) If so, what is the limit? 

(c) Are optimal s-energy configurations uj* n C S d asymptotically uniformly 
distributed on S d 7 

In this paper we show as a corollary to our main results that questions (a) 
and (c) have affirmative answers. Question (b) remains open for d > 2. But 
more interesting is the fact that we can affirm (a) and (c) for a general class 
of (i-dimensional rectifiable manifolds embedded in M. d (cf. Theorem 2.4). 
And for such manifolds, in the case s — d, we give an explicit formula for 
limAr^oo £d(A, N)/N 2 log iV for every d E N. For further background discus- 
sion, see [7,14,15,16]. 



2 Main Results 



In this section we state our main results. Their proofs are given in the sections 
that follow. Let Tid denote rf-dimensional Hausdorff measure in M. d normalized 
so that a d-sided cube with side length 1 has Hd-measure equal to 1. In the 
case d! = d, then H d reduces to Lebesgue measure on R d . 
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Theorem 2.1 Suppose A C W 1 is compact. Then 

S d (A, N) = ^d(B d ) 
NZoN 2 logN H d (A) ' [ ) 

where B d is the closed unit ball in M d . Furthermore, for s > d, the limit 
liniTv^oo £ S (A, N)/N l+S l d exists and is given by 

V ^iM) C '>* (a) 

tfSo N i+ s /d n d {Ay/ d ' 1 ' 

where C s>d is a finite positive constant independent of A. 
Remarks. 

(i) From (9) it is clear that, for s > d, 

£ (U d N) 

C sd = lim Z,' / , (10) 

N^oo ]\[l+s/d ' V ) 

where U d := [0, l] d is the unit cube in ~R d . We further remark that if 
Hd(A) = 0, then the limits in (8) and (9) equal oo. 

(ii) Let B(a, p) denote the closed ball in M. d centered at a with radius p. Then 
the limit with A = B(a, p) in (8) is simply l/p d . 

Theorem 2.2 Let A C W 1 be compact with H d (A) > 0, and oo N = {x ktN }% =1 
be a sequence of asymptotically optimal N -point configurations in A in the 
sense that for some s > d 

E s (u N ) C Std 

w+°/ d ~ n d {Ay/ d ' 1 } 

or 

N^ooN 2 logN H d {A)' K ' 

Let S x denote the unit point mass in the point x. Then in the weak-star topology 
of measures we have 

-E**-* — asN^oo. (13) 



Remark. The convergence assertion (13) is equivalent to each of the following 
assertions: 



(i) For each / continuous on A, 

i N 
im — 

N 



1 N If 

}to^g/(^) = ^4)/ 4 /W^W (14) 
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(ii) For every measurable set B C A whose boundary relative to A has Ttd- 
measure zero, the cardinality \B fl uj n \ satisfies 

\BHu N \ H d {B) 

as N — > oo. (15) 



Theorem 2.3 Let A be a compact set in M, d such that TCd(A) > and A 
{x* N , . . . , x* NN } C A an optimal N point s-energy configuration for A. If 
s > d, there exists a positive constant C = C(A,s,d) such that for every 
N > 2, 

. , * * jc/N'/ d fors>d, 
mml^ x J,N\>{ c/{NlogN y /d fors = ± ^ 

Recall that a mapping (J) : T ^ R d ' , T C R d , is said to be a Lipschitz 
mapping on T if there is some constant L such that 

\ct>{x)-ct>{y)\<L\x-y\ for x,yeT (17) 

and that is said to be a bi-Lipschitz mapping on T (with constant L) 
if 

(l/L)\x-y\ < \(t>{x)-(t>{y)\<L\x-y\ for x,y G T. (18) 
We say that A C W 1 ' is a d-rectifiable manifold if A can be written as 

n 

A=\J MKk) (19) 
fc=i 

where, for each k = 1, . . . , n, Kk C M d is compact and 0^ is bi-Lipschitz on an 
open set G^^> K k . Obviously any compact subset of a rf-rectifiable manifold 
is a d-rectifiable manifold. 

Theorem 2.4 Suppose A C M d ' a d-rectifiable manifold and s > d. If 
s = d, we further suppose that A is a subset of a d-dimensional C l -manifold. 
Then (8) and (9) hold. Furthermore, if Hd(A) > 0, then (13) holds for any 
asymptotically minimal sequence of N point configurations for A satisfying 
(11) or (12). For the case when A is a bi-Lipschitz image of a single compact 
set in M d and Hd(A) > 0, the separation estimates of (16) hold for any optimal 
N -point s-energy configuration. 

Remark. Note that d' does not explicitly appear in (8) and (9) but arises 
only in the norms for the computation of the energy. 

It is shown in [8] that, for the unit interval U l = [0, 1], 

Um £s{u\K) = (s) (s 1} (2Q) 

AT^oo N 1+s 
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where ((s) denotes the classical Riemann zeta function. Hence, using (10), we 
get C S; i = 2£(s) for s > 1. Consequently, Theorem 2.4 gives the following. 

Corollary 2.5 Suppose A is a compact subset of a 1-rectifiable manifold in 
R d ' and s > I. Then 

£ S (AN) _ ^s)_ (21) 



That Corollary 2.5 holds when A is a finite union of rectifiable Jordan arcs 
was shown in [10]. Since a Lipschitz mapping on an interval is absolutely 
continuous, but the converse is not necessarily true, the results in [10] hold in 
cases not covered by Corollary 2.5. On the other hand, Corollary 2.5 applies to 
1-rectifiable manifolds that are not covered by the results in [10] such as, for 
example, when A is the bi-Lipschitz image of a Cantor subset of [0,1] having 
positive measure. 

For the 2-sphere it is shown in [8] that for s > 2, 

l 7™ p ^^-{~^) Cl(s) ' (22) 

where Cl(s) is the zeta function for the hexagonal lattice L consisting of points 
of the form m(l, 0) + n(l/2, y/3/2) for m, n G Z. Consequently (cf. (9)), 

C s , 2 <^j Cl(s). (23) 

It is conjectured in [8] that equality holds in (22) which, if true, would imply 
that equality holds in (23). 

An outline of the remainder of the paper is as follows. In Section 3 we establish 
some basic lemmas on the minimal s-energy of the union of two subsets of lR d . 
Section 4 gives the proof of Theorem 2.1 for the special case when A is the unit 
cube in R d . In Section 5, we verify Theorems 2.1 and 2.2 for almost clopen 
sets in M d . Results on the separation of points in optimal energy configurations 
are established in Section 6. The proofs of Theorems 2.1 and 2.2 for general 
compact sets in IR d is presented in Section 7 and the proof of Theorem 2.4 
appears in Section 8. 



3 Basic Lemmas 



In this section we establish several lemmas that are required for the proofs of 
our main results. First we establish that if A C W 1 is bounded with nonempty 
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interior, then £ S (A, N) grows as iV — > oo with order N 1+s / d for s > d and 
N 2 log JV for s = d. 

Lemma 3.1 Suppose A C M, d is a bounded set with nonempty interior. There 
exist positive constants Co, C\ (depending on A, s, and d, but not on N) such 
that, if s > d, 

CoN i+s/d < N -j < c 1 N 1+s ' d (N > 2) (24) 

and, if s = d, 

C N 2 log N < S d (A,N) < C^ 2 log N (N > 2). (25) 



Proof. We first consider U d = [0, l] d . Let B(x, r) denote the open ball in R d 
with center x and radius r. Then with C = (l/2) d H d (B(0, 1)) we have 

H d (B(x,r) n U d ) > Cr d (26) 

for any x £ U d and r < 1. 

For N > 1, let c^at = {xi, . . . , xjv} be a collection of iV distinct points in U d 
and let 

Ti := min \xi — Xj\. 
Since B(xi, n/2) D B(xj,rj/2) = for 1 < i 7^ j < N, we have 

V /-/AT 

1 = ft d (C/ d ) > '£H d (B(x i , ri /2) n t/ d ) > -]Trf. (27) 
i=i z i=i 

By the Cauchy-Schwarz inequality we have 

/ n \ 2 N N 

n 2 = E^V d/2 <E-fE^ (28) 

\i=l / i=l i=l 

which is known as the harmonic-arithmetic mean inequality. Thus 

JV I N 1 / I \ s l d 

EsM>Y. 7 s-N^ l- d ) (29) 

i=l 1 i i=l JV Vi/ 

/ N 1 1 \ s/d ( N V /d 



N r d I ~ " \ S^ N r d 
=1 iv ' i / \2^j=l ' i , 

where the next to the last inequality follows from Jensen's inequality (or 
Holder's inequality) and the last inequality follows from (28). Since (29) holds 
for any collection uj n of N distinct points in U d , then using (27) we have 

£ s (U d , N) > N 1+s / d C s/d 2- s (N > 2) 
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showing that the lower estimate in (24) holds for A = U d and with C = 

C s/d 2~ s . 

For s = d, the lower estimate in (25) is not so straightforward. For this case 
we shall apply the known result (7). The unit cube U d in F, d can be projected 
onto a subset of S d via the stereographic projection P : R d — > S d defined by 



F{x) = (tx, l-t)e M. a+ \ t = — -. (30) 

Ob | X 

It is easily verified (and well-known in the case d = 2) that for x, y G M, d , we 
have 

|P(x) - F(y)\ = . 2|x ~, y| (31) 



'1 + \x\ 2 <Jl + \y\ 2 
Consequently, for some positive constant C, 

£ d (U d ,N) > C£ d (F(U d ),N) > C£ d (S d ,N), 

and so the desired lower estimate follows from (7). (Later we shall show how 
(7) can be utilized to determine the precise asymptotic behavior of £d(A, N) 
for rf-rectifiable manifolds.) 

For N > 1, let m be the positive integer such that m d < N < (m + l) d . Let 
uon consist of N points selected from U d H (x d /rnj. Then x — y G X d /m for 

u^. For k = (ki, . . . ,kd) G Z d , set HkHoo = max{|fcj|, i = 1, . . . , d} and 
let |k| denote its Euclidean norm. Then 

x^yduj N \ X V\ l<||k||oo<iTi 



k 



< N i+s/d y- y- J_ = jyi+s/d (2j + l) d (2j l) d 



J' TTi 3 s 



j = l \\k\\oo=j J j=i 

<CN i+s/d y 

j=l 3 

For s > d, the sum in the last inequality is bounded from above independently 
of N, while for s — d, it is bounded by a constant times logiV. Thus the 
estimates (24) and (25) hold for A = U d . 

More generally, if A is a bounded set with nonempty interior, then there 
exist r, R > and Xq,Xi G M d such that rU d + x C A C i?£/ d + Xi. Since 
^ s (pf/ d + x, N) = p- s £ s (U d , N) for any p > and x G M d , the estimates (24) 
and (25) follow for A. (For an alternative proof of the upper bound, see the 
proof of Theorem 2.3 in Section 6.) □ 
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Definition 1 Let t s ^ '■ N — > R 6e defined by 



' N l+S ' d ifs>d 
NHogN ifs = d. 



t-a")={"„ i : (32) 



For A C R d ' and a positive integer N, define 

g s4 (A,N) :=£ 8 (A,N)/T 8>d (N) (33) 

and let 

g (A) := liminf G s , d (A,N), g d (A) := lim sup Q„ td (A, N). 

We set 

g 9 , d (A) := lim Q s4 (A, N) 

iv— >oo 

when the limit (as an extended real number) exists. 

If A C R d is bounded and has nonempty interior, then by Lemma 3.1 there 
exist positive constants Co, C\ such that 

C Q <Q s , d {A, N )<Ci 
for N > 2. Hence, g s d (A) and g s d {A) are both positive and finite in this case. 

Lemma 3.2 Suppose s > d and that A and B are bounded sets in R d '. Then 

9 S ,M U B ) > ks,Mr d/S + ljB)- d ' s y S/d . (34) 

Furthermore, ifg gd (A) < oo orifg sd (A) = oo and g sd (B) < oo, and Af is an 

infinite subset ofN and (ujn)n€N is a sequence of sets un C AU B , N G Af, 
such that 

lim = + 9 d (B)-^Y S/ \ (35) 

Af^oo r, h(N) ^- s ' d - s ' d ' 

r '^ nA ' _ 9 -s/ B ^ d/s m) 
N ls/ A ) d/s + 9_J B ) d/s ' 

Remark. If both g (A) < oo and g (B) = oo, then the right-hand sides 
of (34) and (36) are understood to be g sd {A) and 1, respectively; while if 
g (A) = g (B) = oo, then the right-hand side of (34) is understood to be 

— S , & — S , & 

oo. 

Proof. Both g (A) and g (B) are positive since A and i? are bounded. First 
we assume that g (A) and g (B) are finite. Suppose, for JVeN, that cjjv 

— 5 . tit — S . (A, 
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is a set of N distinct points in A U B. Let = un H A and u;^ = un \ ojfj. 
Then 



^(wjv) = E s {ui) + £ s (u;*) +2 



1 



>£ s (c4) + £ s (u;£) (37) 



and hence 



S S (A UB,N)> NA min =N (£ S (A, N A ) + £ S (B, N B )) , 



(38) 



where N A and N B are nonnegative integers. First suppose s > d. Then we 
have 



> liminf min 

Af^oo N A +N B =N 



> liminf min 

~~ N^oo N A +N B =N 



£ S (A,N A ) /N A \ 1+s / d , £ S (B,N B ) / N B \ 1 



l+s/d 



N 



+ 



N 



l+s/d 



N 



-s/d 



N J 



l+s/d 



N . l+s/d 



(39) 



> min 

0<a<l 



9 s /A)a^ + g_jB)(l-ar*/ d [ 



(Note: In the case N A = we set f^)^ = = 0. The case 

N B = is handled similarly.) Let 

F( a ) := £ M (A)a 1+ ^ + gjB) (1 - o:) 1+ ^ (0 < a < 1). (40) 

The reader may verify using elementary calculus that F has a unique minimum 
value F(a*) = (g_ s d (A)~ d / s + g g d (_B)- d / a )~ iA< , where 

This proves (34) when s > d. 

Now suppose (wAr)jvgA/' is a sequence of sets cj^ C iV 6 A/", such that 

(35) holds. We may rewrite (37) in the form 



E B (u N ) ^ E s {u%) (N A \ l+s ' d i E s (co§) (N B 



> 



jyi+s/d — jyi+s/d y jy 



+ 



N i+ s /d y N 



l+s/d 



(N e Af), (41) 



and hence, if j3 is any limit point of the sequence N A /N, N e Af, we get from 
(35) that F(a*) > F(/3). Consequently, (3 = a*, which is equivalent to (36) in 
the case s > d. 

We leave to the reader the remaining cases where at least one of g g (A) and 
g s (B) is infinite. It is helpful to regard separately the cases when N B remains 
bounded or when N B — > oo as A" — > oo. 
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The s = d case of both (34) and (36) follows in a similar manner and is left 
as well for the reader. □ 

For A,B C R d ', let dist(A, B) := ini aeAbeB \a-b\. 



Lemma 3.3 Suppose A and B are bounded sets in M. d ' such that dist(A, B) > 
0. Then 

9sM ^)< {g s , d ( A r d/s + 9sA B r d/s Y s/d • ( 42 ) 

Proof. If ~g s d (A) or g s d (B) equal zero then g s d (A U B) = and so (42) holds 
(note that the right hand side of (42) is understood to be zero in this case). 

Now suppose g s d (A) and g s d (B) axe both positive. Let 5 = dist(A, B) and sup- 
pose N en. Let N A = [a*N] where a* = g s4 ( B ) d/s / {?sA A Y /s + 9sA B ) d/s ) 
and let N B = N - N A . Then 

S S (A UB,N)< S a {A, N A ) + S 8 {B, N B ) + 25~ S N A N B 
< S S (A, N A ) + S S (B, N B ) + 2<T s iV 2 , 



and hence 



G M.iB AT ^ UAN A )r s , d (N A ) Ss(B,N B )T a4 (N B ) 
G s , d (A UB,N)< + - + 25 



2 



t m (7Va) r i)d (N) r s4 (N B ) T s , d (N) T s , d (N) 



Observe that Um Jv ^ 0O r M (^A)/r a>d (iV) = (a*) 1+s / d , lim^oo T Std (N B )/r Std (N) = 
(1 - a*) 1+s / d , and lim^oo N 2 /t s4 (N) = 0. Thus we have 

9sA A us)< <U^)(«*) 1+s/d + M^X 1 - «*) 1+s/d 

□ 



We say that a set A C R d is scalable if A is closed and if for each e > 
there is some bi-Lipschitz mapping h : A — > A° with constant (1 + e) where A° 
denotes the interior of A. For example, a compact, convex set with nonempty 
interior is scalable since, for e > 0, one may choose h(x) = (1 + e) _1 (a; + eu) 
for any fixed u in the interior of A. Similarly, a star-like set is scalable. 



Corollary 3.4 Suppose s > d and A and B are compact subsets of lR d with 
disjoint interiors such that g Std {A) and g s ^ d (B) both exist and A is scalable. 
Then g s ,d(A U B) exists and 

g s ,M ub) = {gsA A r d/s + gsABY d/s y sld ■ ( 43 ) 
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Proof. Let < e < 1. Since A is scalable, there is some bi-Lipschitz mapping 
h with constant (1 + e) such that h(A) C A° and, hence, dist(h(A),B) > 
dist(h(A), A c ) > 0. Then Lemmas 3.2 and 3.3 imply 

(g s ,Mr d/s + 9sAB)- d,s y s ' d < g_JA u b) 

<g s4 (AUB)<g s4 (h(A)UB) 

< {gsM A )T d/s + 9sAB)- d/s Y s/d ■ 

Since g s4 (h(A)) < (1 + e) s g s4 {A), on letting e -> 0, we get (43). □ 



4 The unit cube U d : = [0, l] d . 



In this section we prove that g s ,d(U d ) exists when s > d. We first prove the 
result when s > d by using the self-similarity of U d to obtain estimates relating 
Qs,d{U d , N) at different values of N. For s = d the method is not immediately 
applicable. Instead we use results and techniques developed in [8] and [6] for 
the sphere S d which actually yield gd,d(U d ) explicitly. The proof of the next 
theorem in the case s — d is given separately in Section 4.1. 

Theorem 4.1 For s > d, the limit g s ,d{U d ) := \m\^^ 00 Q s ^{U d , N) exists and 
is finite and positive. Moreover, in the case s = d, 

2ir d / 2 

9d,d(U d ) = n d (B d ) = -— 1 - v (44) 



Proof. Case: s > d. 

We first establish a lemma relating Q St d(U d , N) at different values of N. 

Lemma 4.2 Suppose s > d, 7 e (0, 1) and m is a positive integer. Then there 
is some constant C > (independent of m, N or j) such that 

g s4 (U d , m d N) < 1 - s Gs,d{U d , N) + C(l - 7)-^-^. (45) 



Proof. If m is a positive integer, let I m = {0, . . . , m — \ } d C Z d and, for i e J m , 
set C/ mii := (f/ d + i)/m. Note that ^, d (f/ m ,i, N) = m s Q s4 (U d , N). 

For i G J m , let be a set of TV points in ( / ~fU d + i)/m with minimum energy. 
Let u m d N = Uig/ m Un- If a; G cj^ and y G u;^, then |x - y| > 5 := (1 - i)jm 
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if ||i - j||oo = 1 and \x - y\ > ||i - }\\ 00 /(2m) for ||i - j||oo > 1- Then 

1 



S s (U d ,m d N) < E s {u mdN ) < E S {J N ) + £ 



ie/„ 



\x - y\ ! 



I 



<E 



S s (—U d ,N) + 5~ s 3 d N 2 + 2 s V m s ||i-j||^ 2 
"' lei m 

l!i-j||oc>l 



od j\t2 

ll 1 ~j||oo>l 



m s -f~ s S s (U d ,N) +m' 



V 



< m d+s 7 _s £ s ([/ , N) + 



3 d N 2 



(1-7) 

where K := Skez d \{o} l|k||^f is finite, and so 



+ 2 s KN' 



Q s4 (U d , m d N) < 7- s Gs,d(U d , N) + (l- 1 )- s (3 d + 2 s K) 



— s/nd 



m 



;+d N 2 



T s4 (m d N) ' 



(46) 



Since m s+d N 2 /r S)d (m d N) = N l ~ s / d , the inequality (45) follows from (46) with 
C = 3 d + 2 s K, which completes the proof of Lemma 4.2. □ 

Now suppose e > and < 7 < 1. Let C be the constant in (45) and let iV* 
be such that Q s , d (U d , TV*) < 9 sd (U d ) + 7 s e/2 and C(N*) l ~ s ' d < (1 - 7) s e/2. 
By Lemma 4.2 we then have 



9m := G s ,d(U d , m d N*) < r s {g_ s4 {U d ) + ^e/2) + C(l - ^(N*) 1 -^ 
for any m G N, and hence 

limsup# m < >y- a g (U d ) + e. 



For N > N* , let m N be the greatest integer such that m d N N* < N. Then 

G (U d N) - £s{Ud ' N) < ^(U d ,(m N + l) d N*) _ 

y s A u I* 1 )- Nl+s/d < ( m d N N*y+s/d U + 1 / m N) 9m N +l 

holds for all N > N*, and thus 

limsup£ M (f/ d ,iV) < hmsup(l + l/m N ) s+d g mN+1 < r s g sd {U d ) + e. 

N->oo N^oo ' 

Since this holds for all < 7 < 1 and e > we have that g sd (U d ) < g sd (U d ) 

and, hence, g s ,d(U d ) exists. By Lemma 3.1, g s ,d(U d ) is finite and positive in 
the case s > d, which completes the proof in this case. □ 
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4-1 Proof of Theorem J^.l in the case s = d. 



By Theorem 1.2 we know that g dtd {S d ) exists and is given as in (7). For iVeN, 
let uj* n denote a set of N points in S d minimizing the ci-energy. Also from this 
theorem we have 

KM = HML (47 ) 

n^oo N H d (S d ) { ' 

whenever A C S d is such that the boundary dA (relative to the sphere) has 
Hd(dA) = 0. Such a set A is called an almost clopen subset of S d . 

Lemma 4.3 For N e N, let denote a set of N points in S d minimizing 
the d-energy. If A is an almost clopen subset of S d , then 

= m 

N ^°° T d,d\ N ) n d (8 d ) 



Proof. We first show that for any almost clopen subset K of S d we have 

limsu P I (m < 9dAS d )^jJ-y (49) 

For this purpose we follow the argument given in [8]. Let {x*^}^ denote the 
points of u^f and for each i, set 

U i)N (x) := \ x - x ),N\~ d i xeS d . 

It is shown in inequality (6.6) of [8] that for every r > sufficiently small we 
have 

jj ( * v . g dtd {S d )N\ogN 

Ui, N {x i>N ) < {l _ rdgdd ( S d )) +°r(N) (N-oo). (50) 
Let A(K,N) := {i\x\ N e K} and N K := \A(K,N)\. Then from (50) we get 
E d (u;* N n K) 1 



, N K g d4 (S d ) ( 1 



N (l-r d g dtd (S d )) \logNj- 

Letting N — ■> oo and then r — > in this last inequality we deduce from (47) 
that inequality (49) holds. 
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Now suppose that A C S d is almost clopen (with respect to Tld), and let 
B := S d \ A. For any set K, we put K N := u* N n X. Then, clearly, 

EdK) E d (A N ) E d (B N ) 2 J_ 

WJV) r did (7V) + Td4 (N) + r^(iV) ^ |* - y\*' 1 J 

We claim that, as A" — > oo, the last term in (51) tends to zero. To see this, 
let e > be given and cover dA by an open (relative to S d ) set Q e such that 
7i d (Q e ) < e and 7i d (dVt e ) = (e.g., let Q e be a finite union of open balls). 
Then, since dist(A, B \ Q e ) > and dist(S, A \ Q t ) > 0, it follows that, with 
A t := AnQ t , B e : = Bnfi e , 

2 1 2 1 

limsup — — -r = limsup — — -r (52) 

n^J r d4 (N) x fX N \x-y\ d N ^J r d4 (N) x f- N \x-y\ d 

yeB N ye B N 

< limsup —L-E d (Q»). 

IV—oo Td,d{lV ) 

Since Q e is almost clopen, we get from (49) and (52) that 

,. 2_ 1_ lc , d , H d {tl,) eg i4 (S d ) 

As e > is arbitrary, we have shown that the last term in (51) goes to zero 
as A" — > oo, as claimed. Consequently, 

9dd (S*) = hm ^ = lim f^Q + . (53) 

y ' V ' ^r d ,d(N) ^°°\T dtd (N) T d ,d(N)) 1 ' 

Since A and B are almost clopen and H d (A) + H d (B) = H d (S d ), it follows 
from (49) and (53) that 

and 

E d (B N ) _ d H d (B) 

-oo^ao - 9dAS ] HdWy 

□ 

Lemma 4.4 Suppose A is a compact scalable subset of S d . Then g djd (A) exists 
and 

gdAA)=n d (B d )/H d (A). 

By saying A is a scalable subset of S d , we mean that for every e > there is 
a bi-Lipschitz mapping with constant (1 + e) that maps the closure of A into 
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its interior relative to S d . Clearly the measure of the closure of such a set is 
equal to the measure of its interior, and so any scalable subset of S d is almost 
clopen. 

Proof. Suppose either C = A or C = B := S d \ A. Then (48) holds. We first 
prove that g d4 {C) < H d {B d ) /H d {C). For p > 1 and N e N, let M[N] : = 

lPju$) N \ where l x \ denotes the integer part of x. Let C M[N] := u* M[N] fl C 
and recall that (47) states that \C M ^\/M[N] -> H d (C)/H d (S d ) as N -> oo. 

Then, for iV large enough, we have 



r ft d (S d ) " M[iV] 
from which it follows that there is some N p such that 

N<\C M[N] \ {N>N p ), (54) 
and so £ d (C, N) < E d (cu* M[N] n C) for iV > N p . Thus we have 

- <r\ v Sd{C > N) 

N->oo T did {l\) 

< lim sun UM[iV])^K W nC) 
" r^JV) r<w (M[JV]) 



W d (C) ' 



where (48) was used to obtain the last inequality. Since p > 1 is arbitrary, we 
have g d4 (C) < H d {B d )/H d {C) for either C = A or C = B. 

Next we show g dd (A) > 7i d {B d ) /7i d (A). Let (ajv);veN denote a sequence of 
natural numbers such that liniAr^oo Qd,d(A, ajv) = fi^ d (^)- F° r ^6 N, let 67V = 
\ (T-C d (B) /H d (A))a N ] where \x] denotes the least integer greater than or equal 
to x and let cn = + &at. Since A is scalable, there is a bi-Lipschitz mapping 
h with constant (1 + e) such that h(A) C A . Then 5 := dist(/i(A), B) > 
and, as in the proof of Lemma 3.3, we have 

£ d (S d , c N ) < £ d {h{A) U B, c N ) < £ d (h(A), a N ) + £ d (B, b N ) + 25~ s a N b N 
< (1 + e) d £ d (A, a N ) + £ d (B, b N ) + 25'% 

and thus 

n lad \ ' ( 1 1 \(f/-> / /i \ T d,d\ a N) . n , n , s Td,d{bN) . 25 

Gd,d(S d , c N ) < (1 + e) d Q d4 (A, ajv)— 4 — f + &,d(£, 6jv) + ^—7 



Td,d(c7v) ' ' T^Ctf) log(cTv) 
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Letting N — > oo and then e — > gives 

Using ^ d (S) < H d (B d )/H d (B) and ft d (,S d ) = H d {A) + as well as 

Theorem 1.2, we get fl^ d (^4) > Hd(B d )/Hd(A), which completes the proof of 
Lemma 4.4. □ 

Now we return to the s = d case of the proof of Theorem 4.1. Recall that 
P : M. d — > S d denotes the stereographic projection defined by (30). Let e\ : = 
(1, 0, . . . , 0) G R d and 1 := (1, 1, . . . , 1) e R d . For < 7 < 1, let 

[/ 7 ;= 7 [/<* + ei _ ( 7 / 2 )l = [1 - 7 /2, 1 + 7/2] x [-7/2, 7/2] x • • • x [-7/2, 7/2] 

and let Ay := P([/ 7 ). Note that Ay is scalable, since the mappings a 1— > 
P(rP _1 (a) + (1 — r)ei), for < r < 1, form a family of bi-Lipschitz map- 
pings with constants approaching 1 as r — > 1 that map Ay into its relative 
interior (cf. (31)). Thus (^(Ay) exists and equals Hd(B d ) /Hd(A^) . For x e U 1 
and < 7 < 1/d, we have 1 - 7 < |x| 2 < (1 + 7 /2) 2 + (rf - l)( 7 /2) 2 < 1 + 2 7 . 
Using (31) it follows that for 7 < 1/d, the function h := P -1 is bi-Lipschitz 
on Ay with constant (1 + 7) and such that C/ 7 = h(A 7 ). Then 

9 d , d {U d ) = l d 9d, d (U,) = l d 9 d ,d(h(A,)) < 7 d (l + l) d 9 d ,d(^) (55) 
and, similarly, 

9d/U d ) = > 7 d (l + ir d 9dM,)- (56) 

Since h = P _1 is bi-Lipschitz on Ay with constant (1 + 7), it follows that 
lim 7 ^ + 7" d W <i (A r ) = H d {U d ) = 1 and so 

7<W(A) = l d Hd(B d )/H d (A,) - ^(£ d ) 
as 7 — > 0. Taking 7 — > in (55) and (56) we then have 

W) < ^ d (f/ d ) < ^(f/ rf ) < n d (B d ) 

which completes the proof of Theorem 4.1. 



5 Almost clopen sets in R d 

A Lebesgue measurable set A C M. d is said to be almost clopen (with 
respect to (/-dimensional Lebesgue measure) if 7id{dA) = where dA 
denotes the boundary of A. 
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Theorem 5.1 Suppose A is a bounded almost clopen set in R d . Then g s ^{A) 
exists for s > d and 

g s4 {A) = 9sAU d )n d (Ar s/d - (57) 



Remark. In particular, g s> d(A) = oo if H d (A) = 0. 

Proof. First, if A = ^U d then g s4 {A) = T s g s4 (U d ) = n d (A)- s ' d g s4 (U d ) 
showing that A satisfies (57). Applying Corollary 3.4 inductively, it then fol- 
lows that (57) holds if A is the union of a finite collection of cubes with disjoint 
interiors. 

Next, for n G N, let Q n denote the cubes q in R d with vertices in the lattice 
Z d /n, let A n denote the union of the cubes in Q n that are also contained in 
A and let A n denote the union of the cubes Q n that meet the closure of A. 

Suppose e > 0, then there is an open set V containing dA with Ti. d (V) < e. If 
q G Q n is a subset of A n C\A n , the complement of A n in A n , then q meets dA. 
Since dA is compact and V c is closed, the distance dist (dA, V c ) > 0. 

Let n* be large enough so that diam q < dist (dA, V c ) for q G Q n *. If q G Q n * 
meets dA then q C V and so we have A n n A n C V for n > n*. Hence, 

n d (An) < K d (A n ) < H d (A n ) + e (n > n*) 

showing that 

lim H d (A n ) = lim H d (A n ) = H d (A). 

Since g s4 (A n ) < g s d (A) < g s4 ( A ) < 9s,d(A n ) and (57) holds for A n and A n it 
follows that (57) holds for A. □ 

We say that a sequence ujn C A, N G A/", of sets of points in A is asymptot- 
ically s-energy minimizing on A for s > d if 

Corollary 5.2 Suppose A is a bounded, almost clopen set in M. d , TC d (A) > 0, 
B is an almost clopen subset of A and that ujn C A, N G A/" ; is asymptotically 
s-energy minimizing on A. Then we have 

lim lUN 2 Bl = H d (B)/n d (A). (N G AO (58) 
N nW N 
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Proof. Note that B' := A\B is almost clopen (since dB' C dAUdB). Applying 
Theorem 5.1 to A, B and B' gives 

g.M) = 9sAu d )n d (Ar s/d 

= g s AU d ){H d {B) + n d {B'))- s ' d = (g s , d (B)- d/s + <? M (5T d/s ) ■ 
Then Lemma 3.2 implies (58). □ 



6 Separation 

Proof of Theorem 2.3. For convenience we denote x* N by :rj. For i = 1, . . . , N, 
let 

W-E^K (59) 
Then C/j(xi) < Ui(x) for all i6A Let < 5 < 1 and set 

r := (5H d (A)/(NH d (B(0, l)))) 1 /" 

and 

Dj := B( Xj ,r ), V t := A\\J Dj. 

Then 

Wd(A) > - Nrpi d (B(0, 1)) = W d (^)(l - 5) > (60) 

and we have 

^(xi) < t-^tt / Ui(x)dH d (x) (61) 
= HtW^J^W^ dHd{X) 

- Tup-) § F^F 

Let i? > diam A. It is easy to verify that for < r < 1 and y E A 

I r^—r H dn d (x)<f j\-dn d (u) (62) 

JA\B(y,r) \X - y\ s J B(0,R)\B(0,r) \u\ s 

c s /r s ~ d for s > d, 



< 



c d \og(R/r) for s = d, 



where the positive constants c s , c d are independent of y and r. Using the 
estimates (60) and (62) we get for s > d, 

£ (iV-l) c .(^TO(0,l)))^ £ ^ 

(1 - Oj/l^A) 
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and for s — d, 

U t ( Xi )<k d N log N, (64) 

where the constants k s , k d are independent of iV and %. Finally, since for each 
i = 1,...,N, we have \xi — Xj\~ s < Ui(xi) for i ^ j, inequality (16) follows 
from (63) and (64). □ 

Lemma 6.1 For the closed unit ball B d := 5(0, 1) C M. d , there is a d- 
energy asymptotically optimal sequence (u^Nm of N -point configurations 
ojn = {xi,n, ■ ■ ■ , %n,n} for B d such that for N > 2 

min \x itN - x jtN \ > (2 + \ r d)~ 1 N~ 1/d . 



Proof By Theorems 4.1 and 5.1, we have that g djd (B d ) = 1- For a positive 
integer m let tl m := (±Z) d DB d and for j G (±Z) d let C/ mJ := £[-1/2, l/2] rf +j 
denote the d dimensional cube of side length 1/m (with sides parallel to the 
coordinate axes) and center j. Since 

B(0, (1 - Vd/m)) c (J C/ mJ C B(0, (1 + Vd/m)), 

we have 

(m - Vd) d H d (B d ) < \n m \ < (m + Vd) d H d (B d ). (65) 
Fix k > \fd. If x G C/ m j and |j| > k/m, then |j| > |x| — Vd/(2m) > and so 

0<|j|<2 0<|j|<fc/m fc/m<|j|<2 

" je^zj- 'j' A/m<|x|<2 (|a;| - Vd/(2m))* 1 ' 

0<|j|<fc/m 

l(k-Vd/2)/ 



< 2 d A; d m d + m d / ^ H^AS^dr 

J(k-Vd/2)/m r d 

< 2 d k d m d + m d {\ + \fd/k) d - 1 H d - 1 (S d - 1 ) f 2 -dr 

J(k-Vd/2)/m r 

= 2 d k d m d + m d H d - 1 {S d - 1 ){\ + Vd/kf' 1 log (j^^J^j ■ 
Hence using (65) and the preceding estimate we obtain 

E d (nn = E E iAu* < \n m \ E ,4 

0<ij|<2 

< m d (m + Vd) d H d (B d ) (2 d k d + H d - 1 (S d - 1 )(l + Vd/kf' 1 log ( fc _^ /2 ) 
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Suppose N >2. Now choose m = \(N/H d (B d )) 1/d + Vet] . Then using (65) we 
get 

(m-Vd- l) d H d (B d ) <N < (m - Vd) d H d (B d ) < \n m \. (66) 
Hence we may let u>n consist of N distinct points from Q m . Then 

E d {u N ) < E d (n m ) 



N 2 logN ~ N 2 logN 



< "' 



Xm-Vd-l) 2d J H d (B d )\og((m-Vd-l) d H d (B d )) 
On taking N — > oo (and thus m — > oo) we get 

limsup Ed ^ N "> < Kd-itS*- 1 ) x n, k \d-i _ + 
N^N 2 \ogN- dH d {B d ) ( 1 + Vd W -Cl + Vd//bJ 

for any A; > (here we recall 'H d -\{S d ~ 1 ) = dH d (B d )). Letting k — > oo 
then shows that (cjjv)iveN is (i-energy asymptotically optimal for £> d . Using 
7~Ld(B d ) > 1 and the definition of m, we have m < N l l d (2 + and thus 

min \x-y\ = l/m> (2 + Vd)' 1 N~ 1/d 
which completes the proof. □ 



7 Compact sets 



Proof of Theorem 2.1. Let e > and G be an almost clopen set (since A is 
compact, G could be chosen to be the union of a finite collection of open balls) 
such that G D A and Ti d (G \ A) < e. Then, from Theorem 5.1, 

gJA) > g 8 ,d{G) = gsAU d )n d {G)- s ' d > g s4 (U d )(H d (A) + e)~ s ' d . (67) 

If H d {A) = then (67) shows g sd (A) = g Sjd (A) = oo; if H d (A) > 0, then since 
(67) holds for arbitrary e > 0, we get 

gJA) > g s AU d )H d (A)- s / d . (68) 



We next show g s4 (A) < g s4 {U d )H d {A)- s / d . The case H d {A) = was already 
considered above and so we assume H d (A) > 0. Let 

A* I c Air H d (B(x,r)nA) 

A := {x G A limsup , = , r— = 1}. 

r ^o+ n d {B(x,r)) 
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The Lebesgue Density Theorem (e.g., see [13]) states that H d (A \ A*) = 0. 
For < e < 1, let 

C, := { fl(,, r) | x e A', < r < 1, ^g|g^ >!-<>. (69) 



By the Besicovitch Covering Theorem (cf. [13]), there is a countable collection 
of pairwise disjoint closed balls {Bi := B(xi,Ti)} C C e that covers almost all 
of A* and hence almost all of A. Choose n large enough so that 

(n \ n 

\JADbA = Y / H d (AnB i ) > (l-e)H d (A). (70) 
i=i / i=i 

Let i G {l,...,n} be fixed and let u>n denote an asymptotically minimal 
sequence of configurations for Bi such that 

5 N := min \x-y\> ri(CN)- 1/d 

for some positive constant C independent of %. (Recall Theorem 2.3 states that, 
in the case s > d, any minimal sequence for Bi must satisfy such a separation 
condition while Lemma 6.1 implies the existence of such a sequence in the case 
s = d.) 

For < v < 1/2, let r :— v5n and set 

oj u n := {x G ujn I dist(x, A fl Bi) < r}. 

Then 

B(x,r) n B(y,r) = for x,y £ lon, x ^ y, 

and 

B(x,r)r\A = $ ioi x e u) N \uu u N . 

Since, for any fixed constant less than 1/2, say 1/4, at least this fraction of 
every B(x, r), x G B iy is contained in Bi for iV sufficiently large (and hence r 
sufficiently small), we have 

H d (B t nA c ) >H d { (J B t nB(x,r)) > (l/4)\u N \ u u N \ H d (B(0, l))r d , 
which implies 

\ujn \ oj u n \ < m d {Bi n A c )H d (B(0, l))-\u8 N )- d < ^C—.N, 

v 

where we have used (cf. (69)) U d [Bi n A c ) < eH d [B % ) = erfH d (B(0, 1)). Thus 
\u u N \ = N -\u N \u u N \>N(l-4Cev- d ). (71) 
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If x G uo u N , then there exists y G A fl Bi such that \x — y\ < r. For each x G uo v N , 
let 4>n,v{x) be one such y, and let 

Ajv,„ := {(pN,u( x ) I x e ^n}- 
Note that for x, y G u;^ we have 

|0iv,i/(a;) - <?W(y)| >\x-y\- \4>nA x ) ~ x \~ \0nAv) ~ v\ ( 72 ) 
> (i-2v)\x-y\, 

and so 

E S {\ N ,„) < (l-2u)- s E s (u N ). (73) 
Let M := \N/(1 - ACeu' d )]. Then 

|A M) „| > (1 — ACev~ d )M > N, 

and so we have 

S 9 (A n 5,, iV) < p a (A M ,,) \ / r M (M) \ 



< 



I / r M (M) \ / ^(cjmA 



(l-2i/)- \T s , d (N)J V r M( M )/' 
From the definition of M it follows (even in the case s — d) that 

/ 1 \ l+s/d 

Km t ,am)It,an) = (t^cSr) 

and hence 

^ n BO < (1^7 (ir^)^'*^ (75) 

for any (4Ce) 1 / d < i/ < 1/2. 

Now Theorem 5.1 implies ^(-Bj) = g s ,d{U d )H d {Bi)- s / d and so using Lemma 3.3 
and inequalities (75) and (70) we obtain 

9sM)^9 s , d ({jAnB^ (76) 

1 / 1 \ i+s/d / n \ ~ s / d 



1 

(1 - 2uY 



I / I N 1+s / d 
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for any < e < 1 and any (ACe) 1 ^ < v < 1/2. By first taking e — > and then 
i/^Owe have 

9 s , d (A) < g s , d (U d )H d (Ar s/d (77) 
which combined with (68) completes the proof. □ 

Proof of Theorem 2.2. Let B C A be a measurable set such that Hd(d r B) = 0, 
where <9 r -B := <9.B D (A \ B) is the relative boundary of B. Then A = A\ U A2, 
where A\ :— B U 9 r S, and A 2 := (A \ 5) U <9 r (A \ 5) are compact sets and 
H d (d r (A\B)) = 0. Since H d (A) = H d (A 1 )+H d (A 2 ) we can apply Theorem 2.1 
and Lemma 3.2 to deduce that 

KnA| HJJ3) 

~Jv > HJA) asiV -°°- ( 78 ) 



On writing A = d r B U A \ d r B we similarly have 

\uj N nd r B\ 

— > as N — > 00 

which together with (78) gives (15) and thus (13). □ 



8 (i-Rectifiable manifolds in M. d ' 



In Section 2 we defined the notion of a rf-rectifiable manifold. More generally, 
a set A C M. d ' is said to be a d-dimensional rectifiable set if A is H d 
measurable, Ti d (A) < 00, and TC d -almost all of A is contained in the countable 
union of Lipschitz images of bounded subsets of M. d (see [5], [11], and [13]). 
Clearly, any ci-rectifiable manifold is a <i-dimensional rectifiable set. 

We shall need the following result of Federer concerning ci-dimensional recti- 
fiable sets: 

Lemma 8.1 ([5, 3.2.18], [13, 3.11]) Suppose A C M d ' is a d- dimensional 
rectifiable set and e > 0. Then there exists a countable collection {Ki \ i = 
1,2,.. .} of compact subsets of M. d and bi-Lipschitz mappings tpi : Ki — > M. d , 
i = 1,2,..., with constant (1 + e) such that ^i(Ki), ^2(^2), ^3(^3) •• ■ are 
pairwise disjoint subsets of A that cover Ti d - almost all of A. 

Proposition 8.2 Suppose A C M d ' is a compact, d-dimensional rectifiable set 
and s > d. Then 

9sM) < 9sAU d )n d {A)- s ' d . (79) 

Proof. If Ti d (A) = 0, then the right hand side of (79) is understood to be 00 
and (79) holds trivially. Now suppose < e < H d (A). Let Ki, K 2 , . . . and 
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vpi, ip2, ■ ■ ■ be as in Lemma 8.1. Let n G N be large enough so that 



J2Wi(Ki))>H d (A)-e. (80) 



i=l 



Since ipi(Ki), . . . ,ip n (K n ) are disjoint compact subsets of A, we may use 
Lemma 3.3, Theorem 2.1, (80), and the fact that ifji is bi-Lipschitz with con- 
stant (1 + e) to get 

M^)<2m(0^(W) (81) 

(n \ s/d 

E^«))- d/s J 

/ n \ ~ s / d 

< g s , d (U d )(l + e) 2s \Y,H d {UKi))\ 

<g s4 (U d ){l+6) 2s {H d {A)-tr S/d - 

Since e is arbitrary, (81) shows g 8jd (A) < g s4 (U d )H d (A)- s / d . □ 

Proposition 8.3 Suppose s > d and that A C M. d ' is as in Theorem 2.4 with 
the property that for each e > there is some 5 > such that g sd (B) > 
1/e whenever B is a compact subset of A with TL d (B) < 5. Then g gd (A) > 

g s ,d(u d )H d (A)- s / d . 

Proof. Suppose e > 0. Again let (Ki, ■*/>,), i — 1, 2, . . ., be as in Lemma 8.1. Let 
5 > be such that g sd {B) > (e)~ s l d whenever B is a compact subset of A 
with H d (B) < 5. Let n be large enough so that 

n 

J2Wi(Ki))>Hd(A)-6. (82) 
i=i 

Since A is a Borel set and 1i. d (A) < oo, then 7^^]^ is a Radon measure on R d ' 
(cf. [11, 1.11]). If K C A is compact and e > 0, then there is some relatively 
open set G C A such that K C G and such that 7i d (G) < 7i d (K) + e. 
Furthermore, we may choose G to be Hd-almost clopen relative to A. Indeed, 
if G is not almost clopen then we can construct an almost clopen set Q with 
the same properties as G in the following way. Let C(x,r) = {y E R d ' 
\y - x\ = r}. Since Ti d (A) < oo, the set {r > | H d (C(x,r) n A) > 0} is 
at most countable. Since K G A is compact, there is a relatively open cover 
of K of the form {B(xi, ri) C\ A \ i — 1, . . . , m} where B(xi, r^) n A C G and 
H d (C(xi,ri) n A) = 0. Let (? = U£li r^) H A, then AT C C G, ^ is a 
relatively open subset of A, and H d (Q) = H d (G) < H d (K) + e. 
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Using arguments 2 similar to those in the proof of Theorem 2.1 for s > d (or 
the assumption that A is a subset of d- dimensional C 1 manifold in the case 
s = d), we can find for i = 1, . . . , n a relatively open subset of A such that 
if>i(Ki) C Gi and 

9_s, d (Gi) > fa/MKi))-*'' + e/*Y' /d ■ (83) 

Let G := A \ U™ =1 Then G C A\ U"=i and thus, using (82), we 

obtain Tid(Go) < So and hence 

9 S / G °) > ^ /d - ( 84 ) 
Since ipi is bi-Lipschitz on Ki with constant (1+e) we have, using Theorem 2.1, 

i sd (^m)>a+^r s 9sAKi) (85) 

= (l + e)-'g aid (y d )H d (Ki)-'' d 

> (l + e)- 2s g s4 (U d )H d (UKi)r s/d . 



Since A C Ur=o G ii we again use Lemma 3.2 together with (82)-(85) to obtain 



9_JA)>[T,9jG t )- d/s ) (86) 

Cn n \- s / d 

i=0 i=l / 

/ n \- s / d 

> \2e + (1 + ef d g s4 {U d )- d ' s £ HM(Ki)) ) 

> (2e + (l + e) 2 Vrf(C/ d )- d/ ^(A))~ Vd . 

Taking e — > in (86) then completes the proof. □ 

Proof of Theorem 2.4- Suppose A C R d ' is a <i-rectifiable manifold. Since any 
o?-rectifiable manifold is a <i-dimensional rectifiable set, Proposition 8.2 implies 

g s , d (A) < g s , d (U d )H d (A)- s / d . 

We next show that A also satisfies the hypotheses of Proposition 8.3 which 
will then imply that g s , d {.A) exists and is given by 

9sM) = gsAU d )H d (A)- s / d . (87) 



For further details, see the online addendum arXiv:math-ph/0412053 
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Since A is a rf-rectifiable manifold, we have A = Uk=i 4>k{Kk) where C M. d 
is compact and 0fc is bi-Lipschitz on Kk with constant for = 1, ... ,n. Let 
L := max{L fc | k — 1, . . . , n}. 

Suppose .B is a compact subset of A. For k = 1, . . . , n, let -B^ := B D <pk{Kk)- 
Then by Lemma 3.2 and Theorem 2.1 

> (Es.^)--/') (88) 



>L-g a4 (U d ) 



Vfe=l 



> n-'' d L- 2 'g aid (U d )H d (B) 



-s/d 



from which it follows that A satisfies the hypotheses of Proposition 8.3, thereby 
proving (87). 

Once we have the formula (87), the proof of Theorem 2.2 may be repeated 
without change to show that (13) holds for asymptotically optimal s-energy 
iV-point configurations in A. 

Finally, to prove the separation estimates (16) for an optimal iV-point s-energy 
configuration X* N = {yl N , . . . , y* N N } for A, when A = (f)(K), K C R d , K com- 
pact and bi-Lipschitz on K, we can imitate the argument given in Section 6.1 
for the proof of Theorem 2.3. For this purpose we replace the definition of Ui(x) 
in (59) by 

Ui{x) := §i0(x)-0(x,, w )r 

where x j:N = 4>~ l (y* N ). The details are left to the reader. □ 

Acknowledgements. We thank the referees, S. Borodachov, and J. Brauchart 
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